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ITV ABSTRACT 

on ; 

Q\ . It is shown that one can fit the available data on the ao(980), /o(980), /o(1300) 

and 1Tq(1430) mesons as a distorted ++ qq nonet using very few (5-6) param- 
eters and an improved version of the unitarized quark model. This includes all 
. light two-pseudoscalar thresholds, constraints from Adler zeroes, flavour sym- 

metric couplings, unitarity and physically acceptable analyticity. The parame- 
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ters include a bare uu or dd mass, an over-all coupling constant, a cutoff and a 
strange quark mass of 100 MeV, which is in accord with expectations from the 
quark model. 

> . 

CN . It is found that in particular for the ao(980) and /o(980) the KK component 

in the wave function is large, i.e., for a large fraction of the time the qq state 
t^J- ' is transformed into a virtual KK pair. This KK component, together with a 

similar component of rj'n for the ao(980), and 7777,777/ and n'n' components for 
the /o(980), causes the substantial shift to a lower mass than what is naively 
expected from the qq component alone. 

Qh! Mass, width and mixing parameters, including sheet and pole positions, of 

O ,| the four resonances are given, with a detailed pedagogigal discussion of their 

meaning. 

> 

^ ' 1 Introduction. 

sa : 

As has often been stated in many reviews |lj our present understanding of the light 
meson mass spectrum is in a deplorable state, especially when one considers the vast 
amount of data that has been available already for quite some time, and that QCD in 
principle should solve the hadron spectrum. This is mainly because of the fact that 
the expectations of most "QCD inspired quark models" fail so dramatically for the 
scalar mesons, - the Higgs bosons of the hadronic sector. 

For the other qq nonets such as the 3 Si (p, u, K*, </>), 3 P 2 (a 2 , f 2 , /£), 3 P 1 
(d, fi,KiA, f{), 1 Pi (b\,hi, Kib, h[) and even the 1 So (vr, 77, K, r/) the naive quark 
model works reasonably well as a rough first approximation. Therefore, few authors 
doubt that they should be classified as qq states. Certainly, also here there are some 
"second order effects" such as the observed deviations from ideally mixed states (i.e., 
mass splittings like p — u or A — and mixing angles like <fi — u or — Kib)-, for 
which one has not yet reached consensus as to their full origin. Many authors believe 
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these deviations are mainly due to gluonic intermediate states, while others |2|], includ- 
ing myself, believe the dominant effects come from hadronic loops like — > KK — > uo 
etc. 

But, for the lightest scalars the a (980), / (980), / (1300), and #5(1430) one has 
not even reached a clear consensus as to their true nature. Are some of these qq 
or KK bound states? Or is one of them possibily a glueball? Many authors today 
believe the ao(980) and the /o(980) to be KK bound states]!], [|. This seems, at first, 
to be a natural assumption, since they lie just a little below the KK threshold. If 
so, the 1=1 and the ss state must be sought for at higher masses. And indeed, there 
are now candidates for such states: the a (1450) of Crystal barrel [§]] (or possibly the 
questionable a (1320) of GAMS@) for the 1=1 state, and the LASS / (1525) for 
the ss state. However, flavour symmetric couplings (which works rather well for the 
established nonets) would require that their widths should be at least 500MeV, 
which is much larger than the observed widths of these candidates. Furthermore, 
these can have other interpretations (as radial excitations, meson-meson bound states, 
glueball or threshold effects) like many of the other observed scalars in the much 
too overpopulated 1370-1720MeV region, where many 1=0 candidates 0, /o(1370), 
/o (1450)@, / (1525), / (1590), f 0/2 (1710), do not find a place in the qq quark model. 
One of these could be a glueball, another a deuteronlike state (deuson) etc. In order 
to reach a better conclusion of their true nature it would of course be very helpful if 
at least the lightest scalar qq nonet would be resolved. 

It is fundamental also in many other respects to have a good model of the scalars, 
in particular for understanding chiral symmetry breaking, nuclear forces and of course 
confinement. For chiral symmetry breaking it would be important to know: Where 
is the sigma meson? Is it the / (980), or the / (1300), or must its mass be pushed 



to infinity? How does the Nambu-Jona-Lasinio mechanism] 10] work for the light 



spectrum? Is the pion both the Nambu-Goldstone boson and the 1=1 1 S qq meson, 
as most authors believe, or does one have to look for the qq pion at higher masses, as 
questioned e.g., by Georgi and Manoharfll]. 

For nuclear forces, and for the understanding of the AI = 1/2 enhancement in 
K — > 2tt and K^ decayfl~2l, one would like to have a very light a, in the range of 
600-900 MeV, coupling strongly to tttt; a meson which does not seem to exist. The 
lightest scalar is the /o(980), which behaves more like an ss or a KK state coupling 
weakly to tttt, while the /o(1300) seems too heavy. In this paper I shall suggest a new 
solution to this old question. 

Finally, for confinement the scalars play an important role in building through 
tadpole diagrams the hadronic bag, within which the quarks reside. They are impor- 
tant for the confinement energy, for the qq condensate, and for the difference between 
constituent and chiral quark masses. Thus, they are crucial for the understanding of 
all hadronic masses. 

Many authors have recently studied some of the scalars we discuss here, but gener- 
ally these have tried to fit at the same time only one or two of them and not the whole 
flavour nonet simultaneously. Bugg et al.|13[ and Morgan and Pennington ||14j| have 
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made detailed fits with many parameters to irir amplitudes using general techniques 
of unitarity, analyticity and the K-matrix. Janssen et al. [|l5j studied the a (980) and 



/o(980) in a meson exchange model for irir and irr/ interactions, Achasov|16] studied the 
low energy irir data also with meson exchange and Adler zero constraints. Kaminsky 
et al. [T^J studied the a (980) and / (980) using a coupled channel model conclud- 



ing that the /o(980) is a "KK molecule". Earlier many authors [|n| have discussed 
unitarity and analyticity for the hadron spectrum. 

The results presented here is a very much improved calculation and discussion 
of a short letter |nj 14 years ago. In particular, I have included constraints from 



Adler zeroes, which considerably improve the agreement with data close to thresholds 
without increasing the number of parameters. 

In the following I first discuss the general ingredients and properties of the unita- 
rized quark model (UQM) in Sec. 2. Although this section does contain new results 
and new material to resonance phenomena and the UQM, it can perhaps be skimmed 
by those who just want to understand the fits. It is written in a rather pedagogical 
way and emphasizes some facts, which are often forgotten by many model builders. 
Sec. 3 is the central chapter of this paper, where the actual application and fits to the 
scalar nonet are presented while in the concluding remarks, Sec. 4, some comments on 
these results are discussed. 

2. The unitarized quark model. 

The UQM incorporates unitarity and physically acceptable analyticity with res- 
onances in a way, by which one maintains a simple and transparent physical inter- 
pretation of the introduced resonance parameters. It is a kind of advanced form of 
the classical work of Weisskopf and Wigner pOf . The UQM was applied to many dif- 



ferent hadrons^TJ, and it can explain the signs || and magnitudes of deviations from 
ideal mixing, and many mass splittings. Particularily significant is the large splitting 
between T(6S) and T(45), which cannot be understood in single channel potential 
models (where the predicted splitting is over 50 MeV too small) nor by gluonic ex- 
change. At present there is no other mechanism than hadronic shifts from the loops 
T(nS) — > DD, DD* etc. — ► Y(m5')[]22|, which can account for this large splitting. 

2.1 General formulation 

In the UQM one writes for the partial wave amplitudes (PWA) in "Argand units" 
a factorized matrix form: 

A l3 (s) = T^sXkikj) 1 * =Y, G \ s )^P a p{s)G Pj {s) , (1) 

a/3 

where one sums over the resonance indices a and (3 (for e.g., /o(980), /o(1300) etc.) 
and where % and j denote the two-body thresholds (e.g., irir, KK etc.). The matrices 
G a i{s) include coupling constants (g a i), phase space (and angular momentum) factors, 
and form factors. We write for real s: 
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GUs) = glfi^-F^ S )9( S - s thji ) , s thji = (m Ai + m B f , (2) 

where fcj(s) = [A(s, m\ , m 2 B J /s} 1 ^ 2 /2 is the cm momentum of the two intermediate 
particles Ai and E>i. The form factors F a i(s) are at this stage still quite arbitrary, 
except for the fact that we shall require them to vanish sufficiently fast at oo (In the 
following we generally suppress the index a on the F ai , since in our application in this 
paper we assume this does not depend on the resonance). One expects them to be 
smooth functions of s, which include angular momentum barriers, radial nodes, and 
in principle the left hand cuts. In the quark pair creation model they are given by 



an overlap of the three hadronic wave functions multiplied by a matrix element for 
the qq pair creation. In fact, the Fi(s) include most of the model depencence of our 
scheme. 

For qq resonances the propagator matrix P a p(s) depends on the bare mass param- 
eters mo, a , and on the mass shifts ReH a p(s), and the width-like functions ImU. a p(s), 
which together determine the analytic vacuum polarization functions U Q p(s). 

P^(s) = (mg )Q - s)8 aP + n Q/3 ( S ) , (3) 

The unitarity condition, (A — A') / (2i) = AA^ , takes a very simple form determin- 
ing the imaginary part of H a p(s): 

ImU^(s) = - Y,G m (s)Gl(s) = -Y^g m g^ 1 ^F 2 {s)9{s - s tKi ) . (4) 
i i V s 

Since the functions U a p(s) are analytic functions with only right hand cuts, we can 
write dispersion relations for the real parts ReH a p(s): 



ReU 



7T Js th>1 {S' - S) 



These need no subtractions, since we require that the hadronic form factors Fi(s) 
make ImXl a p(s) go to zero sufficiently fast at infinity because hadrons have finite 
size. Thus the integrals are finite and we need not add any polynomial to ReYi a p{s) 
apart from the a — s term, which we already included in in order to have 
the qq resonances as CDD[Q poles. By defining Rell a p(s) through the dispersion 
relation one automatically satisfies physically correct analytic properties, i.e., one gets 
no spurious poles nor cuts and right asymptotic behaviour. (See discussion in Sec. 2.7. 
below). Thus, once we have a model for the G ai (s) and the bare masses, the PWA 
can be calculated. 

By having the functions IT Q/ g(s) in the inverse propagator one automatically sums 
over all iterated loop diagrams of the Born terms (see the diagrams in Fig. la,b), such 
that the amplitude includes an infinite set of diagrams of the form shown in Fig. lc. 
In Fig. 1 we have included, in addition to the qq resonance terms also contact terms 
(Fig. lb) to be discussed in subsection 2.5. 
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Fig. 1. The Born term for (a) the bare qq resonances, for (b) the contact terms, and (c) the loops 
summed by the functions Yl a p[s) in the inverse propagator. 

2.2 A single resonance 

Some special cases are instructive. For a single resonance P(s) is one-dimensional, 
and assuming real Gi(s) (as can quite generally be done for 2 — ► 2 particle amplitudes) 
the PWA takes the form: 

A = Gi(s)Gj(s) = m B w(^i(s)Tj(s))^ 

13 ml + ReU(s)-s-i'£iG^{s) m 2 (s) - s - im B w^tot{s) 

One recovers a generalization of the familiar Breit-Wigner form in the second expres- 
sion of eq.(||), where one has put I\(s) = Gf (s) /msw- Here m 2 (s) = ml + ReH(s) is 
the "running squared mass", which is given by the bare squared mass plus the gen- 
erally negative mass shift ReU(s). This function is approximately constant only in 
special situations, e.g., if one is far from all thresholds. For S- waves the s-dependence 
is particularily important, since ReH(s) has square root cusps at each threshold, and 
near such a threshold there is a dramatic s-dependence. See Fig. 2a,b where we dis- 
play the running mass and ImU(s) for the Kq and for the ao as obtained in the fit 
presented in the next section. Note in particular the strong cusp at the Krj' threshold. 
Since the A"q(1430) lies sufficiently below this threshold the s-dependence of m 2 (s) is 
here not so crucial as it is for the ao(980) and for the /o(980) resonances, which lie 
essentially right at the KK threshold. One sees from Fig. 2b that for the ao(980) all 



5 



Fig. 2. (a) The real and the imaginary parts of the function m 2 (s) = 171^ + II(s) for the 
-fTg(1430) resonance plotted as a function of y/s. Note the strong cusp at the Krj' threshold and 
that the Krj threshold essentially decouples because of the small coupling constant. The dashed 
curve is S, which crosses the running mass Re[m 2 (s)] at the BW mass. 

(b) As in (a) but for the a (980) resonance. Note that the three thresholds have similar coupling 
strengths and that they lie much closer together than in (a). Therefore the large mass shift for 
ao(980). The dashed curve is S, which crosses Re[m 2 (s)] at the BW mass. 
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three thresholds iri], KK, and nrf lie close to each other. They all contribute to a 
large a (980) mass shift and a strong s-dependence in the running mass. 

2.3 Two and more resonances 

When one considers two resonances the propagator matrix becomes two-dimen- 
sional and has important off-diagonal elements given by Tl a p(s). In general, one 
gets an s-dependent complex mixing angle between the two states, since the mass and 
propagator matrices are diagonalized by a complex orthogonal matrix. This generates 
e.g., OZI violation although the bare states are assumed ideally mixed qq states. 

A special, instructive case is obtained if one has only one threshold, or if the 
coupling constants of the two resonances are proportional to each other, G2%(s) = 
a(s)Gu(s). Such a proportionality can hold when the resonances have the same 
flavour content, say two Kq 's with flavour symmetric couplings. Then both bare states 
have couplings proportional to the same SU^f Clebsch-Gordan coefficients, which are 
given by the general trace formulae: 

g a>l = gTr{M Ai MB i Mc a \c A .c B .c Ca , (7) 

where Ai, Bi, ,C a stand for the three mesons involved at the vertex C a — > A4 + 
Bi, and Ma, is the 3x3 flavour matrix for meson A4. Either the symmetric or the 
antisymmetric trace is to be taken, depending on the sign of the product the three 
charge conjugation quantum numbers. For the ++ — > 0~ + 0~ + vertices, which we 
discuss in this paper it is thus the symmetric trace. Then after some algebra one can 
reduce*the expression ([!]) to: 

A , _ [1 + a\s)]G u .(s)G l3 (s) 



'rriQx — s)(toq2 — s)/{m 2 {s) — s) + II(s) 



where 



a(s) = G 2i {s)/G u {s) = g2iF 2i {s)/[g u F u {s)\ , (9) 
ImU(s) = -[l + a%s)]J2GUs) , (11) 

i 

M (!) = ip r ImWW ■ ("J 

7T Js th S' - S 

(13) 

This shows explicitely how resonances must be "added multiplicatively" because of 
unitarity. Wheras the single resonance form (§) gives one loop in the Argand diagram 

*The result looks almost like magic if one starts from eq. (0), but it is easily found starting with the 
K or even better with the B matrix defined in Sec. 2.7. 



7 



the two resonance form @ gives two loops with a zero in the amplitude at s = rh 2 (s). 
Note that the zero is unshifted, which could be phrased as a theorem: A zero in the 
PWA in the physical region remains a zero after unitarization. 

I find it rather surprising that this quite simple generalization (§) of the BW 
formula to the case of two resonances is not found in the literature, although it could 
be very useful phenomenologically e.g., when studying two resonances in channels 
with nonzero flavour like Kir, where eq. (Q) can hold. 

A special case of eq. @ is also instructive, and this bears some resemblance to the 
"S* effect", i.e. the / (980) discussed below: Imagine two nearly degenerate resonances 
both with large couplings to common channels. Unitarity will shift both masses and 
the phase shift will pass 90° for the first and 270° for the second resonance, close to 
each other in energy. Thus, one gets a narrow resonance structure in the form of a dip 
between two broad bumps. How can this this come about? In fact, when looking at 
the original form, eq. ([]]), the 2x2 mass matrix in the propagator, when diagonalized 
[cf. next subsection 2.4, eqs. (|T7| - |l~9D ] results in a mixing between the two resonances, 
such that one of them nearly decouples from the thresholds, while the other gets very 
large couplings (cf. eq.(^) below). Thus in terms of eigenvalues one has one very 
broad and one very narrow resonance, the latter producing a strong dip in the broad 
resonance bump. Or in other words, two inherently broad resonances can together 
produce one very narrow one! 

For the /o(980) this mechanism is of course not the whole story. Here the couplings 
are not proportional to each other, i.e., they do not satisfy eq.@, and the KK 
threshold plays a crucial role in bringing the two resonances close to each other, 
but loosely speaking a variant of this mechanism is operative. A somewhat similar 
phenomenon appears also for the two axial K\ mesons, which are near 45° mixtures 
of the Kia and Kib, (which belong to the 1 ++ and 1 H nonets) such that one of them 
nearly decouples from K*ir and the other decouples from Kp (cf. Katz and Lipkin|f2|). 
Also in the late sixties there were discussions of similar effects in the connection with 
split resonances [HI (in particular the now well forgotten "split a^). 

For the more general case with two (or more) resonances and many thresholds 
with different kinds of couplings, like (uu + dd)/\^2 and ss states, one does not 
gain much insight by trying to reduce the matrix form of the propagator of eq. (P 
algebraically. The "reduced" formula becomes quite complicated because of the energy 
dependent complex mixing induced. But, fortunately eq. ([[]) is as it stands already 
quite transparent physically, and easy to compute numerically. 

2.4 Breit-Wigner masses and widths, pole positions and resonance mix- 
ings. 

For the single resonance case eq.(^) the Breit-Wigner (90°) mass is given by the 
s- value where m\ + Rell(m s ) — s vanishes or by 

m 2 BW = ml + ReU.(m 2 BW ) , (14) 
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while the widths could be defined as I\ = Gf{m BW ) /rriBw as in eq.(|6]). However, if 
the slope of ReU(s) at the resonance is large, one should correct for this and absorb 
the slope into the s term by dividing both numerator and denominator in eq.(^|) by 
the same term, and define the BW widths by 

r BW _ G}(m BW ) / rriBw v bw _STv BW 

This renormalization of the widths and couplings (which is also familiar in field 
theory) has a clear physical interpretation: Below each threshold the qq pair produces 
virtual pairs of the two mesons, and the probability for such pairs is proportional to 
— j-IIj(s) = — - / ds'ImU(s')/ (s — s') 2 . The qq wave function obtains meson-meson 
components \AiBi >: 

IV> >= [\qq>+Y\-T ReJli ^ AiBi > ]/[l + E(-^n,( S )] . (16) 



i 



Each \A,iBi > component has in configuration space an exponential radial tail, 
whose slope is inversely proportional to how much below the threshold the state is. 
Thus for the ao(980) and /o(980) the size of the KK component grows both in spatial 
size and in absolute magnitude the closer to the threshold from below the resonance is. 
The reduction of the width and coupling in eq. flT5| ) is physically due to the fact that 
only the qq component annihilates directly to tti] respectively tttt. The \KK > part is 
rather inert, since it must first transform into qq near the origin and then into rrr] or 
7T7T. Above the threshold the situation is very different; the KK component vanishes 
since it can simply fall apart, and gives an absorptive part to the wave function. 
The slope of ReU(s) is here generally positive implying that the BW width is in fact 
enhanced. 

In the case when one has many resonances one can define the generalization of 
BW masses and widths by first diagonalizing the mass matrix and propagator by a 
complex orthogonal matrix O(s), which also rotates the couplings which now become 
complex: 



m 



*) = E^Wft^+^Wl^W. ( 17 ) 



2 

diag,a 

a'/3' 



sLM = Do^Wj^, (is) 

a m diag,a\ S ) ~ S V s 

One can then define the BW masses in the multiresonance case as the energies 
where Re[m 2 di aga (s)] — s vanishes. This definition has the advantage that this mass 
in principle is the same in all channels i. We shall here use this definition. Other 
definitions, such as the energies where the phase of each term in the sum flT9] ) is 90°, 



9 



or the energies where their absolute value of each term is maximal, do not have this 
property. But, the analogy with the single resonance case is not simple, because the 
coupling constants g' ai (s) are now complex and energy dependent and furthermore, 
there is a background from other resonance tails. As an example the mixing angle 
between the ss and (uu + dd)/\/2 components in the / (980) or / (1300) is complex 
and is strongly s-dependent. It has a different value at the / (980) than at / (1300), 
and furtermore it will be different when evaluated at the BW mass, than at the pole 
position for the same data and model. 

The pole positions have the advantage that only here does the process factorize 
into production and decay independently of the background. These are determined 
by the complex s value and sheet number, where the the whole inverse propagator 
wIq + n(s) — s vanishes, or more generally where det[P _1 (s)] vanishes. For each 
threshold the number of sheets is doubled, and the sheet number is determined by 
the signs of Im{kj). The same resonance can, in general, have several image poles 
on different sheets, which considerably complicates matters since the same resonance 
can have more than one nearby pole (See Morgan and Pennington [H] and [36 for a 



discussion). In our model for the /o(980) and /o(1300) we have 5 thresholds, which 
means there are 32 sheets, some of which could even have 2 poles each! To find the 
nearest poles one must analytically continue Tl a p(s) defined above at the first sheet 
only just above the real axis. At the first sheet we can calculate U a p(s) by generalizing 
ec[-@ to complex s values on the first sheet (/) by the Cauchy integral around the 
cut on the real x axis: 

K + 1 ,00 IjnU^dx 

7T Js th X - S 

This is discontinuous across the cut. To get to the second sheet (II) one must add a 
term: 

nS(a) = ni p (s) + 2iG al (s)G pi (s) , (21) 

where Gi a (s)Gip(s) is given by essentially the same expression as eq. (£3) or eq.(fj), 
but now without the theta function, and defined also for complex values of s. Near the 
real axis H^q(s) has the opposite sign compared to n£g(s) in the imaginary part above 
threshold. Again below threshold, the cusp in the real part changes sign. Explicitly 
the additional term is given by: 

iG ai (s)Gpi(s) = -g ai gpi^-\(s, rn 2 A .,m 2 Bx )s~ 1 Ff{s) , (22) 
or with the Adler zero constraints and flavour symmetric couplings (See sec. 2.5): 

iG ai (s)Gpi(s) = -7ai7j9i(s - s A ,i)^-X(s,m 2 A .,m 2 Bi ) s _1 i^ 2 (s) , where (23) 
lai = r )\Tr[M.A i M.B i M-c a \+\ ■ (24) 
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In order to get to the third sheet (III) two terms from the first and second threshold 
must be added 

n^/( S ) = Uip(s) + 2iG aX {s)G pi {s) + 2iG a2 {s)G m {s) , (25) 

while to get to the fourth sheet one should only add the second term. In general, with 
the signs of Iraki given by the sheet number (n) one should add 2iG a j(s)G pj(s) for 
each threshold j with negative Imky. 

n^(s) = ny S ) +2iJ2G aj (s)G Pj (s)9(-Imk j ) . (26) 

j 

2.5 Background and t and u channel exchange terms. 

For the case of a nonresonant background term as given e.g., by a contact term 
(Fig. lb) one still has an expression very similar to the one above, but with the m 2 — s 
terms in the propagator replaced by a constant, which could be chosen = 1. But, in 
order to have the same dimensions for the coupling constants of the background as 
for the resonances we put this constant =±M 2 (with dimension GeV 2 ), and with the 
sign allowing for constructive or destructive interference. Then for the contact terms 
we have: 

^contactW = ±M*8 a P + IW(S) • (27) 

The same constraints from unitarity and analyticity still apply as given by eqs. 
(0,fH|), i.e., II Q , / g(s) is given by the same formulas as before although the interpretation 
of the coupling constants now refers to the background. 

We can thus add a background term to each resonance by doubling the dimension 
for the propagator, such that for the case of one resonance and a background one 
writes a 2 x 2 inverse propagator matrix 

The same formulas for n a( g(s) still apply, but since U a p(s) has off diagonal terms 
the background mixes with the resonance in a complex although specified way just as 
in the two resonance case. An important special case appears when the background 
couplings are proportional to the resonance couplings, as is the case when one assumes 
the same flavour structure for resonance and background (g a i,backg = 0-9ai,reson)- In 
this case when one algebraically reduces the dimension of the propagator matrix 
one gets the same resonance formula, but with a modified G a i(s). With the same 
original Fi(s) for resonance and background one needs only to substitute g 2 ai — > 
g\ Jl ± o?{jn\ a — s)/M 2 }. This new linear factor implies that there is a zero in the 
amplitudes, which we identify with the Adler zeroes near s = 0. Thus by adding the 
resonance and contact term with a definite relative weight one introduces the Adler 
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zeroes |27j needed from current algebra. We let the Adler zeroes be at s = SA,i- All 
this is accomplished by the simple substitution: 



9 2 m FKs) - lUs-s A ^F?{s\ or (29) 

k(s) 

G ai (s)GfSi(s) -> j ai jpi(s - s Ati )Fl(s)- ! —^9(s - s th ,i), (30) 

V s 

where we now have introduced dimensionless coupling constants 7^, which are related 
by flavour symmetry eq.(Q). 

One can of course in a similar way include more complicated dynamics coming 
from t- and w-channel exchanges. By making a partial wave projection of the Born 
term one first determines the G Q i(s) and the function, which should replace M 2 or 
Wq — s above, which will contain logarithms etc. Then one sums the iterated higher 
order diagrams by the above unitary formalism. 



2.6 The UQM when two pseudoscalars are produced by some other 
reaction. 

Much of the physics of irn and Kir systems can be learnt from other reactions than 
the two-body reactions discussed above. Such reactions are T' — > Ytytt, K — > -a-atv or 
central production of mr in pp collisions. The UQM can easily be modified in order 
to be applicable also to such processes. All one needs to do (provided the strong 
interactions in the final state is only between the two pseudoscalars) is to replace the 
vertex functions Gpi(s) at one of the two vertices by some other real functions C/3 P (s), 
parametrized in an appropriate way for the production, while the propagator matrix 
including the bare masses and the H a p(s) functions, and Gpi(s) at the second vertex 
for the decay remain as above: 

T ip {s) = G\s) ia P aP {s)C Pp {s) , (31) 

a/3 

This guarantees that the Watson final state theorem, and unitarity in general, 
is automatically satisfied and the Adler zeroes will appear only at one vertex in the 
remaining G a i(s). 

I believe this formalism should be much easier to apply, and physically more 



transparent, than e.g., the formalism of refs. |13|, |14[ where one multiplies the whole 
PWA's by real functions af (s). I have made some initial calculations along these lines, 
and seen that it is easy with appropriate choices of the Cp p (s) to obtain the / (980) 
either as a dip or as a peak in the cross section. But more quantitative applications 
are left for further work. 



12 



2.7 Comparison with the N/D method and the K-matrix. Physically 
acceptable analyticity. 

For a single resonance and one threshold our function G 2 (s) is the N function, while 
our P~ 1 (s) is the D function of the N/D method. For the more general case of many 
channels and many resonances there is no clear connection between the parameters in 
eq. ([!]) and those of a matrix form of N/D[]28|. But, much of the same philosophy of 
N/D methods is also present in the UQM: Given a model for G a i(s) analyticity and 
unitarity constrain the form of the propagator P a p. 

As to the relation between the K-matrix formalism and the UQM, one obtains the 
K-matrix corresponding to the matrix (p]) by simply putting ImU a p(s) = 0: 

KM = { Pl p 3 yi 2 K t] {s) = £ G\s) ia [ml a + ReU a(3 (s) - s^G^s) . (32) 

a/3 

where pi = fcj/ y/s. 

The expression ([!]) is regained from the familiar formula: 

A = K(l - iK)- 1 . (33) 

Note that the mass shift term ReTl a p(s) remains in our K-matrix. In particular, 
thresholds which have not yet opened contribute to ReU a p(s). Thus the UQM can 
be looked upon as a particular way of parametrizing the K-matrix in a way which is 
consistent with analyticity and dispersion relations, whereby one has a simple physical 
interpretation of the parameters. 

By adding the ImH a p(s) terms to the propagator one unitarizes the model and 
sums over the imaginary parts of an infinite series of loop diagrams (Fig. lc). But 
to add the ReU. a p{s) terms is equally important. In many models such terms from 
nearby closed thresholds are omitted, wheras I find it essential to include a complete 
set of nearby flavour related thresholds. Through these one includes the mixing and 
mass shifts of the bare states with the continuum of meson pairs. 

By omitting also the ReU. a p{s) terms from the denominator one obtains the bare 
matrix: 

B*{8) = £ G\s) ta [ml a - s}- l G a3 {s) . (34) 

a 

Since the bare masses and the functions G a i(s) determine the PWA A i3 - through 
the scheme described above, there is a one to one correspondence between the bare 
and the physical masses, once the functions G ai (s) are given. 

One can gain some physical intuition for the UQM formalism, by viewing the 
dispersion relation for ReU(s) as the limiting case of the familiar second order per- 
turbation formula for a mass shift AE = J2i\< at\H\i > \ 2 /{E — Ei). Each piece of 
the continuum shifts the qq state at the same time as the continuum is mixed into 
the state, with an amplitude < a\H\i > /(E — Ej). The sum of the squared mixing 
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amplitudes J2i I < ct\H\i > \ 2 /(E — E/) 2 again corresponds to the above —4-ReTl(s) 
= —- J ImH(s')/[s — s'] 2 ds'. But, in contrast to second order perturbation theory 
the dispersion relation formulas are "exact" in the sense that they solve the coupled 
channel model exactly. 

It is important to emphasize that the function U a p(s) must in general satisfy a 
dispersion relation like the one in eq.(|5|), whereby one automatically has physically 
correct analytic properties for an arbitrary form factor. But in the literature one often 
finds violations of this rule. A few examples should clarify this point: 

If one would put H(s) oc k(s) / i/i, i.e. make it proportional to the relativistic phase 
space factor, (which would be reasonable for Im II (s) if one had pointlike hadrons) 
and use it both for the real and imaginary parts, one would have a spurious pole and 
cut at s = in the physical region. The correct procedure is to calculate the real part 
from a dispersion relation, whereby one gets the Chew-Mandelstam function, which is 
more complicated and has a logaritmic large s behaviour (one needs one subtraction 
constant because of the logaritmic divergence). Then one has no spurious pole nor 
cut at s=0. 

An even simpler example is given by the function (4m 2 — s) 1//2 /s, which also has 
a spurious pole. But, defining the function through the dispersion relation and its 
cut (now one needs no subtraction constant) one gets (4m 2 — s) 1 / 2 / s — 2m/ s i.e., one 
automatically subtracts the spurious pole at s = from the physical sheet. [On the 
second sheet the pole remains, but does no harm: There the function is —(4m 2 — 
s) l l 2 / s — 2m/ s.] As a third example we take a finite cut as given by the function 
[(4m 2 — s)(s — A 2 )] 1 / 2 . The dispersion relation automatically subtracts a polynomial, 
which guarantees that the function vanishes asymptotically in the physical region: 
One gets [(4m 2 - s)(s - A 2 )] 1/2 - s + 2m 2 + A 2 /2. 

These added terms are not small nor insignificant. They alter the model predic- 
tions in a very crucial way, giving e.g., the sharp rise in ReU(s) above the threshold 
(see Fig. 2). But as already mentioned, in the literature one often encounters models, 
which either disregard the real part entirely, or uses some variant of the physically 
unacceptable forms discussed above. Such forms might work within a limited energy 
region, if the spurious singularities are very far away, but would certainly fail if one 
considers several thresholds and a large energy region from threshold to 1.5 GeV, like 
in the present model. 

3. Comparison with data on the lightest scalar mesons. 
3.1 Parameters and form factor 

As discussed in the previous section the PWA's of eq. (P are defined once one 
has a model for the vertex functions G ai (s) in eq. (Q) or ([30]) and the bare masses 
of eq. @. For the latter it is natural to assume an "ideal" and isospin symmetric 
structure such that the bare uu or dd mass mo is a free parameter, while the us and 
ds bare masses are given by mo + m s and the ss bare mass by mo + 2m s , where m s 
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is the bare strange constituent mass. 

For the bare couplings 7 a j we use the OZI rule giving connected flavour symmetric 
couplings of eq. (0). The actual numbers are given in Table 1 for general values 
of the pseudoscalar mixing angle Sp, which measure the deviation from the ideal 
states. The value of 5p is fixed to —54° (which is equivalent to a mixing angle 
Op = 5p + cot -1 y/2 = —19° for the angle measuring the deviation from pure SU3f 
states). Such a value is close to what has been measured in other contexts. (E.g. 
linear mass matrix formulas give 5p = —58°, quadratic —46°, a recent crystal ball 
experiment |Z9"[| quotes —52.2°, Akers et al.|J[J has 5p = —53,7°, while Gilman and 



Kaufman [pl|1 , Baghi et al.[Q and Donoghue et al. find |33| 5p ~ —55°. The results 



are not too sensitive to its actual value as long as it is in this ball park. All the above 
determinations assume the mixing angle to be the same at the 77 as at the 7/ mass. 
This need not exactly be the case, since mixing angles are in general mass dependent; 
cf. our 5s{s) below. This same simplifying assumption for 5p is made also in this 
paper. With our sign conventions the 77 and 77' states are given by: 

\t] > = — sin 5p\uu + dd > /y/2 — cos Sp\ss > , 

|t/ > = + cos 5p\uu + dd > jy/2 — sin5p|ss > . (35) 
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Table 1. The relative Clebsch-Gordan coefficients for the coupling constants as given by 



the trace of eq. (24). The pseudoscalar mixing angle relative to the ideal frame is fixed at 



-54°, while 5/ = cot" 1 y/2 = 35.26° is the ideal mixing angle. 

More generally, one could allow for contributions from disconnected, OZI rule 
violating diagrams for the bare couplings, like 7G = ■y'Tr[A4A i ]Tr[M. B .Ji4c a \ or j' ai = 
•y"Tr[J^4Ai]Tr[M.Bi}Tr[J^Ac a }: which would introduce more parameters, and whereby 
the bare couplings involving the neutral members {v,v' 1 fo? fb) would be different. 
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Fig. 3. The K* branching ratios into Ktt for the resonances on the K * trajectory as a function of 
the squared resonance mass. The branching ratios fall as the square of the function in eq. ( |36| ) with 
/co = 0.63 GeV 2 . A similar exponentially falling behaviour is seen|| for the branching ratios into 
7T7T for the resonances on the fj and pj trajectories. See text. 

Note however, that our physical couplings do not obey the OZI rule, because of the 
nonzero value of the pseudoscalar mixing angle 5p, and because the loops included in 
our formalism generate a scalar meson mixing angle Ss{s). It is possible that most 
of the OZI violation found is of this nature, coming from mixing in the mass matrix; 
therefore for simplicity, we assume as a starting point that the disconnected diagrams 
in the bare couplings can be put = 0. Of course, if e.g. a nearby gluonium state would 
be present, this assumption would have to be relaxed. 

Since the data start only a few 100 MeV above the thresholds, they are not sensitive 
to the exact positions of the Adler zeroes sa, near s = 0, provided these are small 
enough or of the order 0.1 GeV 2 . Thus we put all sa, = 0, except in the case of 
the K7r and the ttti thresholds. For tctt we put it equal to the current algebra value 
sa^ = m "i/2 (although sa^ = would work just as well). For Kn we let the fit find 
the best value, which turns out rather large and negative SA Kn = —0.42 GeV 2 . Current 
algebra including corrections from chiral perturbation theory p4| would predict this to 
be closer to zero. By modifying the flavour symmetry prediction for the K-7T coupling, 
such that it is 10-15% larger than predicted by the Clebsch-Gordan coefficient of a/3/2 
of table 1, one could fit the data with a smaller sa Kv - Thus this determination of sa Kiv 
should not be taken too seriously. 
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Finally for the form factor Fi(s) one assumes the simple Gaussian form: 



F t {s) = exp[-k^s)/(2k 2 )] . (36) 

This, no doubt, is the most drastic assumption of the model. It is used because it 
works, and because forms of this kind are obtained in the quark pair creation model 
(QPCM) P3"| . There, it is given by the overlap of the three qq wave functions of the 
three hadrons at the vertex times a matrix element for the 3 Po quark pair creation. 
With wave functions of Gaussian shape one gets a similar Gaussian factor, multiplied 
by a polynomial in s. The cutoff parameter k is then related to the hadron size R 
through ko = y/E/R, from which with R « 0.7 fm one would estimate ko to be of the 
order 0.6-0.8 GeV/c. In the fit one finds k = 0.56 GeV/c. Including another cutoff 



parameter at the quark pair creation vertex the Orsay group [35) can account for a 
smaller cutoff parameter fc . 

There is another crude phenomenological argument for a form factor of this form: 
If one plots the elastic branching ratios for resonances on the leading trajectory (i.e 
the Kit branching ratios for the Kj, see Fig. 3, or the tttt branching ratio for the 
fj or pj trajectory ) one finds that they fall exponentially with squared mass or 



with k as the square of eq. (pq) with k m 0.63 GeV/c. Of course this is not the 
same thing, but assuming that angular momentum barriers approximately cancel in 
the branching ratio, that the total widths grow much slower, perhaps linearly with 
s, and that each resonance has an exponential form factor similar to ( |36|) in their 
partial widths to exclusive twobody channels etc., one can argue that there is some 
experimental support for such a shape and that the ko parameter is about the same 
for all mesons. 

No doubt, Fi(s) is in reality much more complicated than this. It should in 
principle include all the left hand cuts. But since the left hand singularities lie much 
farther away than the unitarity cuts, which we have included in detail, one can expect 
that the form (RBI) need not be too bad an approximation. 



In table 2 below the parameters of the model are summarized: 

The over all dimensionless coupling constant 7 1.1395 

The uu and dd bare ++ mass m 1.420 GeV 

The constituent bare strange quark mass m s 100 MeV 

The cutoff parameter ko 0.56 GeV/c 

The Adler zero for Kir s Ak -0.46 GeV 2 

The parameter enhancing the 777/ couplings (3 1.6 



Table 2. The parameters of the model. In addition, the pseudoscalar mixing angle 5p 
is fixed to a conventional value of —54°, and the Adler zero for 7T7T is given the current 
algebra value while for the remaining channels the Adler zeroes s Ai are put 

= 0. Other constants used in the model include the pseudoscalar masses and the 
Clebsch-Gordan coefficients of Table 1. 
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Fig. 4. (a) The Ktt S-wave phase shift (solid curve) and absorption parameter Tj * 100 (dotted 
curve) as functions of the Kn center of mass energy when the model is fitted to the LASS data[^6| 
shown by the round dots. The error bars of the data are of the same order as the dots in the figure. 
This fit fixes four of the six parameters of table 2: 7, TTLq + m s , ko, and Sa Kw - 
(b) As in (a) but for the absolute value of the S-wave amplitude. 



18 



Fig. 5. The Ktt Argand diagram for the same fit as in Fig. 4. The numbers indicate the Ktt 
center of mass energy in MeV. 



3.2. The K*(U30) and the Kn S-wave. 

It is natural to start the data comparison with the Ktt S-wave, since this is the 
simplest to understand, having only one resonance, the Kq{1430), and in addition the 



experimental data are rather good. The comparison of the fit to the LASS datap6 



is shown in Figs. 4a,b. The older data of Estabrooks et al.|37j are very similar but 
have larger error bars. The error bars of the LASS data are of the same magnitude as 
the dots in the figures. Fig. 4a shows the S-wave phase shift, while Fig. 4b shows the 
absolute magnitude of the amplitude (in Argand units). The corresponding Argand 
diagram is shown in Fig. 5. As can be seen the model has no difficulty in fitting 
the data rather well. These data essentially fix four of the parameters in Table 2: 
7, m +m s , ko and SA K7r - We say essentially, since some of the parameters are strongly 
correlated, in particular ko and 7, and in practice because of the time consuming 
numerical integration in eq. (||) one must keep ko fixed. 

The parameters of the ^(1430) also turn out to be rather close to the conventional 
ones (cf. Table 3). The BW mass is at 1349 MeV where the phase shift passes 90°, 
and the BW width is 398 MeV with negligible correction from eq. (|1^), since ReU(s) 
is almost flat at 1350 MeV. The nearest pole is also where one normally expects it on 
the third sheet at [Re(s)] 1 ^ 2 = 1440 MeV, with an imaginary part —Im(s)/m = 320 
MeV (Since the Krj coupling almost vanishes this threshold is very weak and, in fact, 
there is an image pole on the second sheet at almost the same position). 
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Fig. 6. (a) The 717] S-wave Argand diagram as predicted by the model when four of the parameters of 
table 2 are fixed by the Kit data of Fig. 4 and m s = 100 MeV. Note the sharply falling amplitude 
after the KK threshold due to the strong absorption to this channel, which partly explains the 
narrow do (980) shape. The numbers are the center of mass energies in MeV. 

(b) As in (a), but for the phase shift (solid curve) and absorption parameter (dotted curve) plotted 
as a function of 7TT] center of mass energy. 
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Fig. 7. The a (980) resonance shape (solid curve) as predicted (not fitted) by the same model and 
parameters as in Figs. 4 and 6, and compared with the data of Armstrong et al.j38| from central 
production of TTTTTj in 300 GeV pp collisions. The background (dashed curve) is a rough estimate 
of the reflection of other contributions to the experimental distribution and is added to the model 
prediction. Note that in spite of the inherently large width as measured by — IrnHys) /m ~ 400 
MeV or by the imaginary part of the pole position the predicted peak width is rather narrow ~ 100 
MeV. 
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3.3. The a (980) and the nr/ S-wave. 



When turning to the irr] data for the ao(980) there is only one new parameter, m s , 
whose value one expects to be close to 100 MeV. And in fact, the fit allows for this 
value, therefore it is fixed at 100 MeV. The value 100 MeV is slightly lower than what 
is usually quoted in the light sector (e.g., the mass spittings of — K* or K* — (p, uj) 
is about 120 MeV.) But, part of this splitting is renormalized by the unitarity effects 
we include, and in the fits to the other light multiplets^I] values close to 100 MeV 



were found. A better comparison is with the D s — D splitting of 99.1 ±0.6 MeV, since 
here unitarity effects are smaller. 

The data for the ao(980) is much poorer than for the ^(1430). In Fig. 6a the 
Argand diagram for the irr] — > in] S-wave is shown, and Fig. 6b shows the phase 
shift with the absorption parameter. Fig. 7 shows the the prediction for the ao(980) 



peak compared with the data of Armstrong et al.[38] of centrally produced r]7rn in 
pp collisions. It is significant that the model predicts a rather narrow peak at the 
right mass. This low mass comes out right since the mass shift is the largest at the 
threshold, and the three thesholds for the ao: nr], KK, ttt]' (Fig. 2b) lie close to each 
others, while the corresponding ones for K$: Kir, Kr], Kr]' , (Fig. 2a) are much more 
spread out in energy and the middle one, Kr], almost decouples. Again the rather 
narrow peak structure of about 100 MeV seems at first very surprising, considering 
the very large imaginary part (cf. in Fig. 2b, Imll(s)/m »s 400 MeV). This is due to 
several contributing effects, which we list in their order of importance: 

i) The large negative slope of ReU(s) below the KK threshold, which according 
to eq.(|6|) renormalizes the BW width by a large factor, which is close to 5 at ^/s = .98 
GeV; 

ii) The strong absorption above the KK threshold (cf. the Argand diagram of 
Fig. 6a) making the amplitude drop fast in magnitude above this threshold; 

iii) The Adler zero at s = reduces the amplitude somewhat near the threshold. 

See also the discussion of pole position, Sec. 3.5, for another way of understanding 
the narrow peak. 



3.4. The / (980) and / (1300) and the tttt S-wave. 

Finally we turn to the much more complicated tctt S-wave shown in Fig. 8 together 
with the CERN-Munich [[39|1 , Cason et al.|4(|, Grayer et al. ff4TT| and some lower energy 
experiments^]. Here we have two resonances with complex energy dependent mixing. 
Now essentially all of our parameters except one are already fixed by the Kit and tit] 
data. But even without new parameters, the prediction for the phase shift is not too 
bad below 1.1 GeV. The additional parameter (3 is needed in order to fit the phase 
shift above this energy, which otherwise comes out about 40° too low. We chose this 
to be a factor in front of the rjr]' couplings, which thus violate the flavour symmetry 



22 



Fig. 8. (a) The Argand diagram of the 7T7T S-wave as predicted (not fitted) by the model when 4 of 
the 6 parameters of table 2 are fixed by the Ktt data of Fig. 4 and m s = 100 MeV (as in Figs. 6-7) 
and /3 = 1.6. The numbers are the center of mass energies in MeV. 

(b) The same prediction as in (a) for the 7T7T phase shift (solid curve) and absorption parameter 
?7 * 100 (dashed) as function of the center of mass energy, and compared with the CERN-Munichp9[ , 
Cason et al|I(J, Grayer et al.[^4| and some low energy experiments |42| . 
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predictions of Table 1, being increased by 60%. This determines our final parameter 
P = 1.6. Of course this does not necessarily mean that this is the right choice. Equally 
well one could think of other effects, left out of the model, which could be faked by 
/3. These include higher thresholds which enter at about the same energy as rji]' in 
the 1.5 GeV region such as vector- vector or pseudoscalar-axial thresholds etc., or an 
effect due to the higher mass scalar resonances seen in this region. 

In Figs. 9 and 10 the running masses are shown before, respectively after diagonal- 
ization of the mass matrix. Note that in the region of the two resonances the real parts 
of the two diagonal running masses are nearly degenerate (Fig. 9a) due to the mass 
shifts from the KK threshold. Therefore the mixing angle (Fig. 10c) between the 
two resonances changes rapidly at the KK threshold from being nearly ideal 5s ~ 0° 
below 900MeV to 5s ~ —35.26° where the states are close to being SU3f eigenstates. 
Note that the large value of —Im(mf) results in a very large BW width (880 MeV) 
of the first BW resonance ("the a"). Fig. 10a,b shows the eigenvalues, with the BW 
masses given by the crossing points of the real parts with s. In Fig. 10 the / (980) 
does not appear as a BW resonance, but it is hidden under the complex and rapid 
energy dependent mixing. It is another linear combination of the two states, which 
appears when one analytically continues to the nearby poles. 

3.5. Pole positions and BW resonance masses and widths. 

When one asks the question, where are the poles of the resonances in this model 
with large imaginary parts, one is tempted to say: Who cares where the poles end 
up in the complex multisheeted structure of the s-plane! Each resonance may require 
many poles on different sheets for a full description of its pole structure. On the 
other hand one has a very simple bare spectrum, and a rather simple description of 
the data. The bare masses are strongly shifted nonperturbatively to some s values 
and sheet numbers, whose positions depend sensitively on the precise positions of 
where the thresholds are. Do we really learn anything fundamental from knowing the 
precise pole positions, apart from having numbers to put into the PDG tables? I am 
inclined to answer: not very much, they are of secondary importance! Of course, the 
experimental data, resonance shapes and threshold positions are very important in 
order to find the correct solution with a reasonable bare spectrum together with the 
other model parameters. But the actual positions of the poles depend more on the 
positions of the thresholds and on Clebsch Gordan coefficients, than on the model 
parameters. 

The actual pole positions do tell us something of the analytic structure of the 
solution found, and can throw some light on resonance shapes. E.g., the fact that 
the nearest ao(980) and /o(980) poles are found on the second sheet, although their 
Re (sp i e ) are above the KK threshold, throws some light on why these are narrow. 
Normally with two open thresholds the nearest resonance pole is expected on the 
third sheet. The fact that it often turns up on the second sheet explains partly their 
shapes: The nearest place to the pole in the physical region is at the threshold, since 
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Fig. 9. The real (a) and imaginary (b) parts of the functions — m oa^a/3 + n a ^(s) for 

the (uu + and the SS resonances before diagonalization of the mass matrix. Note that the 

SS running mass Re m| 2 (s) dips strongly towards the (uu + dd) / \/2 running mass Re Tn\-y(s) 
near the KK treshold. 25 



Fig. 10. The real (a) and imaginary (b) parts of the eigenvalues of the mass matrix m\,(s) after 
diagonalization, and (c) the scalar mixing angle 5s(s). The crossing points with S (dashed) in (a) 
give the BW masses. Note that 8s{s) is almost real at the energies of the BW masses, and has 
a large imaginary part in the KK threshold region. For low energies the mixing is nearly ideal 
(8s = 0°), wheras above 1.1 GeV one has nearly pure SU3f eigenstates (8s = —35.26°). At the 

lighter BW mass y/s = 0.86 GeV one has an extremely broad (880 MeV) near (uu + dd)/y/2 
BW resonance (the "<7") while at 1186 MeV one has a near SU3f octet BW resonance. On the 
other hand, when one analytically continues to the nearest pole, these combine linearily to give one 
very "narrow" almost pure SS state [the /q(980)], which lies on the second (not third) sheet. See 
text for discussion. o<? 



one must go around the KK cut to see the pole. Therefore, a narrow peak appears at 
the threshold in spite of a very large imaginary part of the pole position, and in spite 
of the fact that Re(s poie ) is about 90 MeV above the KK threshold. Of course, as 
explained above one can also understand the narrow a (980) by the arguments given 
earlier by the structure of the amplitudes in the physical region at Re(s) + it. 

In Table 3 I list the Breit-Wigner masses, widths (in MeV) and mixing angles 
together with the nearest pole positions, the Riemann sheet of the poles and mixing 
angles at the poles. Note that the a and f lie on the second sheet although the KK 
threshold is open at the pole mass, as was already discussed above. Similarily the 
second fo lies on the 3rd sheet although the 7777 threshold is open, and one normally 
would expect it at the 5th sheet. The a width in the T BW column is the peak width 
since it is not possible to define a BW width. 

The mixing angle 5 S is such that at the BW mass of the first / , whose BW 
width is enormous, 880 MeV, 5$ is small, i.e. one has an almost pure [uu + dd)/^/2 
state with a little mixture of ss. The sign of this mixing such that it is mainly an 
SU3f singlet state. This BW resonance we interpret as the "a". On the other hand, 
when one analytically continues to the pole with lowest mass one has another linear 
combination of the two resonances: A narrow almost pure ss state (but with a large 
ilm5s, P oie = ^39°). This is our / (980). See also our discussion in Sec 2.3 of how two 
broad resonances can make one narrow. The heavier /o(1300) is close to an SU3/ 
octet both at the BW mass and at the pole position. Also the mass of the heavier 
state is similar both at the pole (1202 MeV) and when measured as a BW resonance 
(1186 MeV). 
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II 
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Table 3. Breit-Wigner masses, widths (in MeV) and mixing angles together with the 
nearest pole position parameters, Riemann sheet number and mixing angles at the poles. 
Except for the Ss,bw at the light fo the three other mixing angles are given with respect to 
the ss state. The heavy fo is thus close to an SU3f octet, for which 5s = —35.26°, while 
the lighter BW fo is a near (uu + dd)/\^2 state (the a). But at the poles one has other 
linear combinations of the two resonances such that the lighter is a narrow almost pure ss 
state, the / (980). 

The ao width in the Tbw column is the peak width as it is here not possible to define a 
BW width. See text. 

Our solution suggests a novel way to resolve the old question: Where is the a 
meson? From the low energy side the nn interaction looks very much as if there 
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would only be an extremely broad (880 MeV) BW resonance with a mass of 860 
MeV, and with a composition near (uu + dd) / a/2. The fact that the running mass is 
not constant makes the effective a mass slightly vary depending on what energy one 
is sensitive to. It is interesting that this prediction of the a — > tttt width is close to the 
values predicted in the linear a model and in extensions of the Nambu-Jona-Lasinio 
model|l. 

On the other hand as one approaches the KK theshold the two resonances expe- 
rience a dramatic s-dependent mixing, such that at the KK threshold the relevant 
linear combination appears as an almost pure ss pole [the / (980)]. We find the pole 
to be a little above the threshold and on the second sheet, with its precise position 
rather sensitive to model parameters. But, it is always a few MeV from the KK 
threshold. 

4. Concluding remarks. 

In 1972 Feynmanp| wrote in connection with the problem of how to treat virtual 
hadronic loops, needed by self consistency in the quark model, that "it is the most 
important problem in the theory of strong interactions", but continued that "no cal- 
culation of such virtual strong interactions in any problem has ever been successful" . 
He thought this "dead end is a result of lack of imagination of how to get further". 

No doubt, these were strong words, and certainly we have made a lot of progress to 
understand virtual hadronic states since 1972, but the progress has been rather slow, 
possibly because too many theorists really thought that one had reached a "dead 
end" , and looked for other approaches guided by QCD. Thereby the nonperturbative 
methods and related problems of the sixties were "swept under the rug". 

Today it seems more and more obvious that attempts where one, starting directly 
from the QCD Lagrangian, tries to solve the nonperturbative light hadron spectrum, 
has come to another "dead end", in spite of the heroic attempts of e.g., lattice QCD. 
This seems all the more obvious, with the present solution of the scalar nonet at hand, 
which I believe to be close to the true solution, although certainly improvements can 
and should be done. To disentangle all the complicated threshold singularities, shifting 
the poles to some strange corners of the the many-sheeted structure of the complex 
s plane, and then try to understand why some peaks like the a (980) or / (980) are 
narrow in spite of their large couplings, seems to be an impossible task to resolve 
directly from the QCD Lagrangean. It would require inclusion of fermionic loops and 
too much computing power, if one would not know how to proceed through some 
intermediate steps. 

Such an intermediate step which I believe should be useful, is to use QCD to 
calculate, or at least understand qualitatively, the parameters used in the UQM, i.e., 
the overall coupling constant 7, the cutoff parameter k and the bare masses, when 
the mass shifts and couplings to the lightest thresholds have been unfolded. 

A few comments are in order as to the meaning of the bare masses found in this 
paper. The value of ttiq depends on the cutoff parameter fc and on the number of 
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thresholds, which have been included. The larger the cutoff and the more thresholds 
included in the UQM, the larger will tuq be. Furthermore, if one also includes tadpole 
graphs with hadron loops, these have no imaginary part, but would contribute to i?ell, 
and shift m down, partly cancelling the mass shift from the thresholds. Thus the bare 
mass m is model dependent. On the other hand, the bare constituent strange quark 
mass m s is less dependent on distant singularities, since in the mass difference m^ss) — 
mo(uu) = 2m s the above effects essentially cancel. Similarily OZI rule violating 
mixing between resonances (8s) is less dependent on distant thresholds, since when 
one sums over complete sets of F-coupled and D-coupled flavour related thresholds 
their contribution to off diagonal elements in the mass matrix is small [[21] , [45]. Only 
the overall bare mass scale mo can depend strongly on distant singularities. 

The model we have found for the scalar nonet is, in principle, not very different 
than what one would assume for any of the other qq nonets with finite decay widths. 
For these one must always have an input mass scale mo and a strange quark mass m s 
and an overall coupling parameter which, together with flavour symmetry describe the 
finite widths. Also the cutoff parameter must, in principle, exist by self consistency, 
although one implicitly assumes that things do not depend on it. Our solution has 
the same four parameters plus only two additional phenomenological parameters, (/3 
and SA Kir )i which may be just artifacts of the things left out of the model. Apart from 
the P parameter flavour symmetry is broken in our scheme most importantly by the 
pseudoscalar mass differences and by the small m s . 

An important feature of the model is that the ao(980) and the /o(980) have very 
large components of KK (and also some ryV or rji] etc.) in their wave functions. 
We estimate the KK component to be dominant in the ao(980) near the peak being 
about 4-5 times larger than the qq component. This does of course not mean that 
the a (980) or the / (980) are KK bound states; they owe their existence to the qq 
component. But, it is important to note that this will make these states look "as if 
they were KK states". E.g., their widths to 77 will be reduced by this same factor 
of 4-5, and thus a reduced 77 width does not prove that they are KK molecules, as 
has beed argued by Close and Barnes and collaborators [Q. One should not compare 
the KK molecule model with a too naive qq model, without the large virtual KK 
component in the wave function. 

Finally a natural question many readers certainly ask: Why is it that the scalar 
mesons are so much more sensitive to nonperturbative unitarity and analyticity effects 
compared to the other light nonets? And a related question: Why has not the solution 
presented here not been found previously? 

To answer the first question one observes two special features for the scalars: 

(i) The over all squared coupling coupling constant 7 2 is about 5-6 times larger 
for the scalars than for the vector nonet. This can be seen already by comparing 
the ^(1430) width of 287 MeV, which is 5.7 times larger than the K* width 
of 50 MeV. A better comparison taking into account phase space etc. is found by 
determining a corresponding dimensionless 7 2 parameter for the vectors from the 
width expression: (mT)^* = 3/27 2 (4fc 2 )/c/ v /i, which also gives a factor of about 5 
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smaller r j 2 . Theoretically this is roughly what one also expects from spin counting or 
SU6jy symmetry, which would predict a factor of 3. Thus the pseudoscalar thresholds 
are about five times stronger for the scalar mesons than for the vector mesons and 
consequently unitarity effects are much more important. 

(ii) The second feature which makes the thresholds more important for the scalars 
is that these are S-wave thresholds. Therefore one has strong square root cusps in 
n(s) with a discontinuous first derivative. For P-wave and higher angular momenta 
this is smoothened out by the angular momentum factor k 2e . Thereby the scalars are 
much more sensitive to the threshold positions, especially when the resonance turns 
up near a threshold as is the case for the a (980) and for the / (980). 

As to the second question, why this solution has not been found previously, one 
notes that (disregarding our first short note[^9|), no one has tried to fit simultaneously 
the whole nonet taking into account all the light pseudoscalar thresholds, putting in 
physically acceptable analyticity properties etc. 
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Figure captions 



Fig. 1 The Born term for (a) the bare qq resonances, for (b) the contact terms, 
and (c) the loops summed by the functions U a p(s) in the inverse propagator. 

Fig. 2. (a) The real and the imaginary parts of the function m 2 (s) = m 2 ^ + IT(s) 
for the (1430) resonance plotted as a function of yfs. Note the strong cusp at the 
Kr)' threshold and that the Kr\ threshold essentially decouples because of the small 
coupling constant. The dashed curve is s, which crosses the running mass Re[m 2 (s)] 
at the BW mass. 

(b) As in (a) but for the ao(980) resonance. Note that the three thresholds have 
similar coupling strengths and that they lie much closer together than in (a). Therefore 
the large mass shift for ao(980). The dashed curve is s, which crosses Re[m 2 (s)] at 
the BW mass. 

Fig. 3. The K* branching ratios into Kir for the resonances on the K* trajectory 
as a function of the squared resonance mass. The branching ratios fall as the square 
of the function in eq. ( |36|) with k = 0.63 GeV 2 . A similar exponentially falling 
behaviour is seenJJ for the branching ratios into tttt for the resonances on the fj and 
pj trajectories. See text. 

Fig. 4. (a) The Ktt S-wave phase shift (solid curve) and absorption parameter 
i] * 100 (dotted curve) as functions of the Kn center of mass energy when the model is 
fitted to the LASS data|3B| shown by the round dots. The error bars of the data are 



of the same order as the dots in the figure. This fit fixes four of the six parameters of 
table 2: 7, m + m s , k , and sa Kiv - 

(b) As in (a) but for the absolute value of the S-wave amplitude. 

Fig. 5. The Kit Argand diagram for the same fit as in Fig. 4. The numbers 
indicate the Kir center of mass energy in MeV. 

Fig. 6. (a) The ttt] S-wave Argand diagram as predicted by the model when four of 
the parameters of table 2 are fixed by the Ktt data of Fig. 4 and m s = 100 MeV. Note 
the sharply falling amplitude after the KK threshold due to the strong absorption to 
this channel, which partly explains the narrow ao(980) shape. The numbers are the 
center of mass energies in MeV. 

(b) As in (a), but for the phase shift (solid curve) and absorption parameter (dotted 
curve) plotted as a function of nr) center of mass energy. 

Fig. 7. The a (980) resonance shape (solid curve) as predicted (not fitted) by 
the same model and parameters as in Figs. 4 and 6, and compared with the data of 
Armstrong et al.[|38| from central production of tttti] in 300 GeV pp collisions. The 
background (dashed curve) is a rough estimate of the reflection of other contributions 
to the experimental distribution and is added to the model prediction. Note that in 
spite of the inherently large width as measured by —ImH{s)/m ~ 400 MeV or by the 
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imaginary part of the pole position the predicted peak width is rather narrow 100 
MeV. 



Fig. 8. The Argand diagram of the nn S-wave as predicted (not fitted) by the 
model when 4 of the 6 parameters of table 2 are fixed by the Kir data of Fig. 4 and 
m s = 100 MeV (as in Figs. 6-7) and f3 = 1.6. The numbers are the center of mass 
energies in MeV. 

(b) The same prediction as in (a) for the inr phase shift (solid curve) and absorption 
parameter rj * 100 (dashed) as function of the center of mass energy, and compared 
with the CERN- Munich [|31S|] , Cason et al|40||, Grayer et al.[H] and some low energy 
experiments [^] . 

Fig. 9. The real (a) and imaginary (b) parts of the functions m 2 a Js) = m^ a 5 a ^ + 
n Q/ g(s) for the (uu + dd)/\/2 and the ss resonances before diagonalization of the 
mass matrix. Note that the ss running mass Re rn\ 2 {s) dips strongly towards the 
(uu + running mass Re mf x (s) near the KK treshold. 

Fig. 10. The real (a) and imaginary (b) parts of the eigenvalues of the mass matrix 
m 2 a ,(s) after diagonalization, and (c) the scalar mixing angle 5s(s). The crossing 
points with s (dashed) in (a) give the BW masses. Note that 5s(s) is almost real at 
the energies of the BW masses, and has a large imaginary part in the KK threshold 
region. For low energies the mixing is nearly ideal (5s = 0°), wheras above 1.1 GeV one 
has nearly pure SU3f eigenstates (5s = —35.26°). At the lighter BW mass yfs = 0.86 
GeV one has an extremely broad (880 MeV) near (uu + dd) / ^/2 BW resonance (the 
"cr") while at 1186 MeV one has a near 5773/ octet BW resonance. On the other 
hand, when one analytically continues to the nearest pole, these combine linearily to 
give one very "narrow" almost pure ss state [the /o(980)], which lies on the second 
(not third) sheet. See text for discussion. 
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